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ABSTRACT 

We discuss the tight interconnection between microlensing optical depths, flattening 
of dark haloes and low-to-intermediate redshift baryonic census. By analysing plots 
of the microlensing optical depth as a function of galactic coordinates for different 
values of axis ratio q of the galactic MACHO halo, we have shown that observations 
are best described by a flattened halo with q ~ 0.6. There is no dynamical obstacle 
for such a choice of global halo shape. Both extremely flattened q ~ 0.2 and spherical 
q ~ 1 haloes have several difficulties, although not of equal severity. Consequences 
of such flattening for the cosmological density fraction contained in MACHOs are 
considered and comparison with mass in low and intermediate-redshift Lya forest and 
other plausible reservoirs of gas is discussed in context of a unified description of the 
evolution of baryonic content of the universe. 

Key words: Galaxy: halo - Galaxy: gravitational lensing - dark matter - galaxies: 
evolution 



1 INTRODUCTION 

The question of the nature and properties of the baryonic 
dark matter (BDM) is one of the most active fields of re- 
cent astrophysical research (Hegyi & Olive 1986; Persic & 
Salucci 1992, 1998; Ashman 1992; Gnedin & Ostriker 1992; 
Richstone et al. 1992; Carr 1994 and references therein; 
Wasserman & Salpeter 1994; Flynn, Gould & Bahcall 1996; 
Graff & Freese 1996), because of its impact on the var- 
ious branches of the modern astrophysics and cosmology. 
Microlensing searches have proved to be one of the most im- 
portant tools for investigation of properties of the halo of the 
Milky Way (Paczyhski 1986; Griest et al. 1991; De Riijula, 
Jetzer & Masso 1992; Paczyhski et al. 1994; Aubourg et al. 
1993; Sackett & Gould 1993; Gould 1994, 1996; Alcock et al. 
1996, 1997a, b, c; Ansari et al. 1996). Comparison of theoret- 
ical models and microlensing (ML) data has already yielded 
intriguing results and insights (Gates, Gyuk & Turner 1995; 
Steigman & Tkachev 1998). Under the Copernican assump- 
tion that the Milky Way is a typical zero-redshift L, galaxy, 
it is natural to ask what consequences recently discovered 
Massive Compact Halo Objects (MACHOs) have for the 
global picture of evolution of the baryonic structure in the 
universe. 

By MACHOs we denote present-day collapsed objects 
residing in the halo of the Milky Way (and, by a Copernican 
assumption, haloes of other normal spirals) that are bary- 
onic in nature and/or origin. Thus, we exclude hypothetical 
primordial mini black holes (e.g. Canuto 1978) and exotic 



non-baryonic aggregates (e.g. Kolb & Tkachev 1995; Eich- 
ler 1996). The motivation behind this is the crucial impor- 
tance of the Big Bang nucleosynthesis (BBNS) constraints 
for cosmological density fraction Qb = Pb/PctH in baryons 
(Yang et al. 1984; Walker et al. 1991; but see also Gnedin 
& Ostriker 1992; Hata et al. 1996). Objects detectable only 
through microlensing searches may present an important, 
and indeed a dominant, item in the total baryonic census. 

We hereby intend to extend the discussion of an im- 
portant recent paper by Fields, Freese & Graff (1998; here- 
after FFG). While they approach the problem of cosmolog- 
ical density fraction S^macho contained in collapsed halo 
objects from the point of view of mass-to-light ratio and 
MACHOs as products of stellar evolution, we take a slightly 
different approach to achieve similar ends. In our view, mi- 
crolensing observations give the most significant key to the 
spatial distribution of baryonic matter in haloes of large spi- 
ral galaxies like the Milky Way. Under the assumption that 
our galaxy is typical for the normal galactic population at 
zero redshift, we apply integration of spherical and various 
flattened model MACHO distributions over the luminosity 
function in order to determine cosmological contribution of 
these dark haloes. We also extend the discussion of FFG on 
the comparison between MACHO and Lya forest mass den- 
sity, especially in light of recent observational indications 
that most of the low-redshift Lya forest is associated with 
galaxies (Spinrad et al. 1993; Lanzetta et al. 1995; Chen et 
al. 1998; Yahata et al. 1998). This conjecture is crucial to 
the understanding of transition between different types of 
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baryonic dark matter, which does not seem to occur later 
than the differentiation of galactic structure itself. This is 
in accordance with the cooling flow-type models of galaxy 
formation (Nulsen & Fabian 1995, 1997), and offer simple 
interpretation of the best available cosmic baryon census 
(Fukugita, Hogan & Peebles 1998, hereafter FHP), strength- 
ening the assumption of FFG that MACHOs were formed 
before the Lya systems. Thus, our paper is complementary 
to FFG in several ways. 

Global flattening of the galactic dark halo has been re- 
cently considered in detail by a number of authors (Binney, 
May & Ostriker 1987; Griest et al. 1991; Sackett & Gould 
1993; Frieman & Scoccimarro 1994; Nakamura, Kan-ya & 
Nishi 1996). Wc use and update their results, especially ex- 
cellent discussions of Sackett & Gould (1993) and Frieman & 
Scoccimarro (1994) which, unfortunately, were written be- 
fore the bulk of the existing microlensing data came, in order 
to emphasize the tight relation between the global shape of 
the baryonic halo and its total mass, with all its cosmolog- 
ical implications. The goal of this paper, accordingly, is (1) 
to update the discussion of discrimination of various halo 
models according to microlensing optical depths produced, 
and (2) to investigate the influence of realistic global flat- 
tening, constrained in such a way, on estimates of baryonic 
cosmological density and consequences for the evolution of 
the baryonic content of the universe. Thus, inclusion of more 
phenomena, like the global flattening, present another step 
towards a unified picture of the physical processes of rele- 
vance to the baryonic matter. 



2 BARYONIC DARK MATTER IN THE 
GALAXY: A SHORT OVERVIEW 

From shape of the Milky Way rotation curve (RC) (Mer- 
rifield 1992) one can see that a huge amount of mass still 
has to be identified. The difficulties in the determination 
of the RC led to uncertainties in the most important pa- 
rameters such as the galactic constant Rn, which represents 
the distance to the Galactic centre, and the circulax speed 
at the Solar radius, vq (Merrifield 1992; OlUng & Merri- 
field 1998; Sackett 1997). Although the lAU 1986 standard 
values are _Ro = 8.5 kpc and uo = 220 kms^'^ some recent 
estimates allow the smaller values: i?o = 7.1 ± 0.4 kpc and 
Wo = 184±8kms"^ (Oiling & Merrifield 1998). In this paper 
we adopt the value Ro = 8.5 ± 0.5 kpc (Feast & Whitelock 
1997) based upon an analysis of HIPPARCOS proper motion 
of 220 Galactic Cepheids and t;o — 210 ± 25kms"'^ that in- 
cludes the best values from the Hi analysis (vq = 185 kms~^) 
and the estimated value based on the Sgr A* proper motion 
vo = 235 kms"^ (Sackett 1997). 

Without going into the discussions about the content 
of the dark matter in the halo, we only state here that one 
part (presumably smaller) has to be in the baryonic form. 
Namely, cosmic nucleosynthesis predicts that (Turner 1996; 
FHP): 

0.0062 < flB < 0.026 (1) 

where Qb is the universal baryonic mass-density parame- 
ter {Qb = ps/pcrit = 8TvGpB/SHi) and 0.4< h <1.0. The 
"silent" h is used in parametrization of the Hubble constant 
Ho = 100 h kms~^ Mpc~^. Recent measurements of the pri- 



mordial deuterium abundance (Buries & Tytler 1998) give: 

Q.Bh^ = 0.0193. (2) 

Using the simplest dynamical estimate of the mass of 
the Galaxy (Kepler's third law) one can obtain (Roulet & 
MoUerach 1997, and references therein): 

iWdyn ^ ^ =^ 

- 5-6 X 10" ( ^^J"" , \ ( r^^^ Mo, 
\^ 220 km/s y \^ 50 kpcy 

where M(r) is the mass interior to Tmax, Wc is the measured 
rotational velocity and r is the radius within which most of 
the dynamical mass of the Galaxy is located. For the entire 
visible matter (stars, interstellar and intraclustcr gas) one 
can obtain cosmological density fraction (Persic & Salucci 
1992; see also Bristow & Phillipps 1994): 

fivis w 2.2 X 10"^ + 6.1 X 10"* h~^-^. (4) 

The discrepancy between the Eqs. (1) and (4) represents the 
so-called problem of missing baryons. Various types of such 
dark baryonic material have been suggested: gaseous clouds 
of plasma or neutral atoms and molecules, snowballs or icy 
bodies similar to comets, stars, planets, white dwarfs, neu- 
tron stars and stellar or primordial black holes (e.g. Hegyi 
& Olive 1986; Peebles 1993). 

On the other hand, the dynamical mass in the Eq. (3) 
corresponds to the cosmological density parameter (if the 
Milky Way is a typical L* galaxy) . 

(This result depends on details of the Galaxy luminosity 
function which will be discussed in the Section 5). The dis- 
crepancy between the Eqs. (3.1) and (4) can be regarded 
as a formulation of the general dark matter problem on the 
galactic scales. The necessity for the dark matter is empha- 
sized by severe limits on mass-to-light ratio in the Local 
Group imposed by deep blank sky surveys (Richstone et al. 
1992; Hu et al. 1994; Flynn, Gould & BahcaU 1996), as weU 
as with huge dynamical mass for the Milky Way inferred by 
Kulessa & Lynden-Bell (1992) and Lee et al. (1993). 

The mass in the halo is dominated by the matter that 
is not, at least easily, detectable. According to the recent 
observations of satellite galaxies (Zaritsky et al. 1997) dark 
haloes extend to much larger radius than 50 kpc. So, one 
can formally write: 

f^dyn ~ 0.1 ^15nvis. (5) 

3 MICROLENSING AND FLATTENING 

In searches for the BDM content the method of microlens- 
ing has so far proved successful. Its name derives from the 
fact that lensing of distant objects is made by bodies with 
masses characteristic of a star or planet. Although the theo- 
retical development of this idea started in 1964 (e.g. Peebles 
1993, and references therein), it was the seminal paper by 
Paczyriski (1986) that showed that one can search for ML 
events in the Milky Way halo if it is made of stars or brown 
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dwarfs. Rapid development of observational and computer 
technology led to the detection of a significant number of 
ML events (e.g. Mellier, Bernardeau & Van Waerbeke 1998). 
Searches have been directed towards Large and Small Mag- 
ellanic Clouds (LMC and SMC) (Alcock et al. 1996, 1997b; 
Ansari et al. 1996; Palanque-Delabrouillc et al. 1998), Galac- 
tic bulge (Kiraga & Paczyhski 1994) and M31 (Crotts 1992, 
Crotts & Tomaney 1996; see also Gould 1994). 

All these surveys give results concerning two important 
parameters: masses of the intervening lenses and the number 
of lenses within the Einstein radius around the line of sight 
to a lensed source — the optical depth. In Table 1 we give 
the targets observed, names of the appropriate survey, mass 
ranges of the lenses, and corresponding optical depth. 

It is known from the work of Sackett & Gould (1993) 
that instead of the equation for the mass density in a spher- 
ical halo: 



(6) 



(where r is the Galactocentric radius, Voa is the asymptotic 
circular speed of the halo, Vc is the core radius of the halo 
and Rt is the truncation radius), one should use the general 
formula for flattened halo: 



p(r) = 



1 



tanip vie 



0{Rt - 



(7) 



where 



'' + tan^i/) {z denotes height above the Galac- 



tic plane). Here the flattening parameter is introduced: 
cosi/' = 1 = c/a, i.e. its cosine determines the shape of the 
halo En. The En notation is related to g as g = 1 — n/10. 
Following Sackett & Gould (1993) we write the following ex- 
pression for the estimate of the optical depth as a function 
of Galactic coordinates / (longitude) and 6 (latitude): 



D 



{D-L)LdL 



{2Ro cos I cos b)L + {l + sin^ b tan^ ijj)L^ 



(8) 

whore we use Ro — 8.5 kpc and Tc = 5 kpc (e.g. Alcock et 
al. 1996). Now we integrate this equation and take D — 50 
kpc (for LMC), D = 63 kpc (for SMC) and D = 770 kpc 
for M31 (Binney & Tremeine 1987). While Sackett & Gould 
(1993) use values for q starting with q = 0.4 (shape E6), we 
will start with admittedly extreme value q = 0.2 (shape E8) 
suggested by some theories, like the halo molecular clouds 
or the decaying daxk matter (see the discussion below). 

There are several other lines of reasoning suggesting a 
high degree of halo flattening in spiral galaxies. One is for 
long time suspected (e.g. Ninkovic 1985, 1991; Bahcall 1986; 
Binney et al. 1987) flattening of the Population II subsys- 
tem, which may be a consequence of the residual rotation, or 
more probably, global flattening of the gravitational poten- 
tial created by dark matter. A detailed discussion of these 
questions was given by Binney et al. (1987), who attributed 
the flattening to a highly anisotropic velocity dispersion ten- 
sor, under the assumption that the velocity ellipsoid of halo 



objects near the Sun has the same shape as that of the ex- 
treme Population II subsystem. Their calculation show that 
the parameter q for the halo isodensity contours has to be 
in the range q — 0.3 — 0.6. Wysc & Gilmorc (1989) review 
different arguments based on analyses of halo star counts 
for flattened halo and suggest values q = 0.6 as the opti- 
mal one. The same is, with additional data, repeated in the 
extensive review of Gilmore, Wyse & Kuijken (1989) where 
the authors concluded that star count data deflnitely favor 
0.6 < q < 0.8 models. One should note that widely used 
Bahcall-Soneira model of the Galaxy indicates q = 0.8, al- 
though it is realized that it may be too conservative, even 
for stellar subsystem (Bahcall 1986). These results supersede 
older arguments for spherical stellar component (e.g. Oort & 
Plant 1975). As far as other galaxies are concerned, observa- 
tional data are still extremely scarce, but it is indicative that 
very recently, the observations of the gravitational lens sys- 
tem B1600-I-434, consisting of two spiral galaxies (Gl and 
G2), where G2 is a barred one, suggest that it has g^O.4 
(Koopmans, de Bruyn & Jackson 1998). The flattening of 
the M31 halo was discussed and appropriate corrections to 
the mass estimates were considered by Ninkovic & Petro- 
vskaya (1992). 

As pointed out by Sciama (1990), a further theoreti- 
cal virtue of the halo flattening idea is connected with the 
Oort limit. Since the amount of the dark matter per unit 
surface area of the disc is determined by the rotational ve- 
locity of the disc, any flattening of the vertical dark mat- 
ter distribution must be compensated for by an increase in 
the density of daxk matter in and near the galactic plane 
(i.e. in the disc). Therefore, the amount of the dark matter 
near the Solar system, traditionally associated with the Oort 
limit is reduced, and could be even brought down to zero 
(see also Binney et al. 1987). This is appealing, since many 
recent results show incompatibility with the large quanti- 
ties of local unseen matter (e.g. Binney & Merrifield 1998). 
Two different approaches yielded different values for the lo- 
cal density, po near the Sun (the Oort limit): Kuijken & 
Gilmore (1989) found that po = 0.10 M© pc"^, while Bah- 
call (1984) found po = 0.21 M© pc~^. However, it seems 
that recent HIPPARCOS measurements favour lower value 
(po = 0.11 ± 0.01 Me pc"^') (e.g. Kovalevsky 1998). This 
value is not in disagreement with the assumption of maxi- 
mal disc (Sackett 1997). 

It is, also, necessary to consider the behaviour of gas 
distributed in the halo. If the seminal idea of Bahcall & 
Spitzer (1969) of extended gaseous haloes of normal galax- 
ies producing narrow absorption features in the spectra of 
background objects is correct, as indicated by recent low- 
redshift measurements (Bergeron & Boisse 1991; Spinrad et 
al. 1993; Lanzetta et al. 1995, Chen et al. 1998), then the 
distribution of gas could tell us something about the shape 
of the gravitational potential. It is not a simple problem at 
aU (e.g. Barcons & Fabian 1987; Pitts & Tayler 1997), but 
some results are quite suggestive. In an important recent pa- 
per. Ranch & Haehnelt (1995, hereafter RH95) have shown 
that for the most plausible values of Lya cloud parameters, 
the conclusion that their axial ratio (thickness/transverse 
length) is less than 0.25 is inescapable. This conclusion does 
not depend on the exact choice of model for Lya clouds, and, 
if the observations quoted above are correctly interpreted, 
would mean that the gaseous haloes are also flattened by the 
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Table 1. Targets in diflferent ML surveys, the mass ranges of the lenses and optical depths. 
References: (1) Alcock et al. 1997b, (2) Renault et al. 1997, (3) Alcock et al. 1997a, (4) Crotts 
& Tomaney 1996, (5) Palanque-Delabrouille et al. 1998. 



Target 
LMC/SMCi 

LMC^ 
Gal. bulge'' 
M31'» 

LMC/SMC^ 



Survey 
MACHO 

EROS 

MACHO;DUO;OGLE 

KPNO 

EROS2 



Mass range 
0.3 - 0.5 Mq 



0.08 



0.1 Mq 
0.6 Mq 
1.0 Mq 



2.6 



+8.2 



Mr. 



Optical depth 



: 2.9^ 



X 10- 



TLMC - ^-^-0.9 
TSMC = 1.5 - 3 X 10 



7 
7 

X lO-'^ 
X 10"® 



7"LMC = 0.82+j!,;g 

_ Q n+1.8 
'bul — 1.2 

TM31 = 5 - 10 X 10-« 
TSMC = 3.3 X 10-7 



same amount. One should keep in mind that such absorp- 
tion studies probe only "a tip of an iceberg", since these 
objects are ionized to extremely high degree, and may as 
well contain dominant part of the baryonic density in the 
Eq. (1). We shall return to discussion of these questions in 
the Section 6. 



4 OPTICAL DEPTHS FOR FLATTENED 
HALOES 

Bearing this in mind, we solve the integral in the Eq. (8) and 
give estimate for r in several cases of particular interest: 

• Optical depth r(Z,6) in the parametric space, with the 
parameter q fixed in steps of 0.2, i.e. q = 0.2, q — 0.4, 
q = 0.6, q = 0.8 and g « 1. 

• Optical depth r for different targets: LMC, SMC, M31 
and Galactic bulge (bar) in order to see what value of q 
determines the optical depth that is closest to observed 
value in the corresponding survey. 

We hereby present several such three-dimensional plots. In 
Figs. 1, 2 and 3, values of the optical depth r are plotted 
against galactic coordinates I and b for three chosen halo 
models, with q = 0.2, 0.6 and 1 (spherical halo), respec- 
tively. These axial ratios were chosen as representatives of 
extremely flattened (like Pfenniger et al. [1994] or the DDM 
models discussed below), moderately flattened (like those 
favored by Bahcall-Soneira models or other dynamic discus- 
sions) and unflattened models (usual approximation). One 
can use such plots (on the same or smaller angular scales) 
in order to choose the observing direction where the optical 
depth reaches the maximal value. Such an example is shown 
in Figs. 4 through 6, where we have plotted optical depth 
as function of coordinates I — 280. °5 and b — — 32.°9 of the 
Large Magellanic Cloud. 

After solving the integral in the Eq. (8) for given values 
of the parameter q we looked for the values that match the 
optical depth obtained in various surveys. We note that, in 
general, the best agreement can be attained with 0.2^ q 
^0.6. Specifically: 

1. For the case of the LMC, that has been studied rather 
well, the measured value of the optical depth based 
upon the sample of 8 events is r = 2.9lg g x 10"^ (Al- 
cock et al. 1997b) while we find that for q = 0.5 we have 
Tx3x W'^ (see Fig. 5). 

2. For the case of the SMC, that is studied less thoroughly, 
the optical depth is estimated as r = 1.5 — 3 x 10"^ 
(Alcock et al. 1997c). Our results show that the model 
in the Eq. (8) gives the value r^4 x 10-^ for q « 0.5 
and above. 




Figure 1. Optical depth for sources located at Z3 = 50 kpc in a 
flattened model with q = 0.2. 




Figure 2. The same as in Fig. 1, except that model with inter- 
mediate flattening, q = 0.6 is chosen. 

3. For the case of the galaxy M31 we found r ~ 5 x 10"® 
which is an accordance with the estimates 5 — 10 x 10"® 
(Crotts & Tomaney 1996), under the assumption that 
g>0.2. 

4. Determining r towards the Galactic centre is more com- 
plicated and we will not discuss it here. We only state 
that using the model in the Eq. (8) we can estimate 
the halo contribution to the ML rate towards Galac- 
tic centre which is between r « 5 x IQ-^ (g — 0.2) 
and r ~ 1.6 x 10-^ {q ~ 1); the estimated range 
for the total optical depth towards Galactic centre is 
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Figure 3. The same as in Figs. 1 and 2, but for a spherical (q = 1) 
case. 

(r = 3.9li:^ X 10"^) (Alcock et al. 1997a). This dis- 
crepancy might be the consequence of the fact that the 
Galactic disc is maximal: if this is so, then one can ex- 
pect a higher optical depth towards the bulge due to disc 
lenses and lower values of optical depth towards targets 
in the halo, that is measured, cf. Table 1 (Sackett 1997). 
This nice illustration of the possible discrimination be- 
tween various halo shapes, is presented in Figs. 7, 8 and 
9 (the same cases as above). 
From these estimates and studies summarized in Table 1, it 
seems that both the spherical {q ~ 1) and extremely flat- 
tened [q ~ 0.2) dark halo can be ruled out with high sig- 
nificance. The optimal value for q is the intermediate one 
q ~ 0.6. Recent research shows that it is not uncommon 
case with spiral galaxies, not only polar-ring ones (Sackett & 
Sparke 1990; Sackett et al. 1994; Oiling 1995). This has sev- 
eral far-reaching consequences, like implausibility of molec- 
ular halo or DDM models, which, in the simplest forms, pre- 
dict strong flattening of dynamically important halo compo- 
nent. 

We also present ratios of optical depth toward LMC and 
SMC as functions of the flattening parameter q in Table 2. 
This is significant, since variation in tlmc/tsmc is claimed 
to sufficiently clearly discriminate between various models of 
flattening (Sackett & Gould 1993). Our results are in agree- 
ment with the more detailed discussion of this problem by 
Frieman & Scoccimaro (1994). 

The necessity of having many more lines of sight for 
microlensing survey besides LMC, SMC, Galactic bulge and 
Andromeda galaxy, led several research groups to consider 
globular clusters as ML targets (Rhoads & Malhotra 1998, 
Jezter et al. 1998, Gyuk & Holder 1998). The original idea 
came from Paczynski (1994) who proposed 47 Tuc (mon- 
itored by OGLE collaboration) and M22. Gyuk & Holder 
(1998) and Jezter et al. (1998) composed lists of appropriate 
clusters for microlensing survey, with corresponding optical 
depths, stating that by using globular clusters it would be 
possible to distinguish the flattened models of the galactic 
halo much easier, and it would allow better determination of 
halo structure parameters, such as the power-low index, or 
the core radius. Special advantage of that method is avoiding 
self-lensing events present in a SMC survey, for example, and 
fluctuations in halo density due to clumps of stars in tidal 




Figure 4. Optical depth toward LMC for very flattened (g = 0.2) 
case. 




Figure 5. The same as in Fig. 4, except that q = 0.6 case is 
shown. 



tails and intervening dwarf galaxies like recently discovered 
Sagittarius galaxy. 

We briefly note that very recently the EROS collabo- 
ration (Derue et al. 1999) found that in the four directions 
towards the Galactic spiral arms the average value for the 
optical depth is equal r = O.SS^Q jg x 10~^. If one compares 
values obtained for their four targets and the predicted val- 
ues of r from Fig. 2, one could again conclude that the halo 
is moderately flattened, i.e., q ~ 0.5 — 0.6. 

With globular clusters distributed throughout the Mil- 
ky Way's halo, this method will further improve on the num- 
ber of ML events and make it possible to constrain Milky 
Way's dark halo shape and content more precisely. 

It is of foremost importance to continue theoretical 
studies of these optical depths, even if we are not in po- 
sition to find other possible sources of background light for 
microlensing experiments. Further improvement in number 
of events will enable precise measuring of abovementioned 
optical depths, and using their ratios to constrain possible 
halo shapes. 
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Table 2. tsmc/'''lmc ratio for different 
values of the parameter q. Optical depths, 
T, are expressed in units of 10"'^. 



9 


(tlmc) 


(tsmc) 




0.2 
0.6 
1.0 


3.8554 
4.0733 
4.2184 


5.2082 
5.6973 
5.6670 


1.3509 
1.3987 
1.3434 




Figure 6. The same as in previous two Figures, but for spherical 
(q = 1) model. 





Figure 8. The same as in Fig. 7, for moderately flattened halo 
(g = 0.6). 




Figure 7. The halo contribution to the microlcnsing optical 
depth toward the Galactic Centre for q = 0.2. 



Figure 9. The same as in previous two Figures, for the spherical 
case. 



5 TOTAL MASS OF MACHOS AND 
COSMOLOGICAL DENSITY PARAMETER 

The total mass of a MACHO halo of an L, galaxy (as typi- 
fied by the Milky Way) in a model characterized by the Eq. 
(7) is 



M{q,RT) = 3.648 x 10" 



q arccos q 
r dzdr 



(9) 



Mr. 



where all lengths are in cm, and has a fixed value. We 
shall briefiy discuss on the variation of core radius below. 

The total mass of an L, galaxy is shown in Fig. 10 
for q varying between 0.2 and 1, and a truncation radius 
Rt between 40 and 60 kpc. Horizontal plane represents the 
dynamical value of the total mass inferred from satellite 
studies within much larger radius of ~ 230 kpc by Kulessa 
& Lynden-Bell (1992). The choice of interval for the trun- 
cation radius is relevant not only because it incorporates 
the "canonical" value of 50 kpc for the size of MACHO 
haloes (FFG; Alcock et al. 1996, 1997a), which is reasonable 
from the point of view of empirical detection of microlcnsing 
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Figure 10. The mass of the Galaxy in units of IO^^Mq as a func- 
tion of the truncation radius, Rj- and the flattening parameter, 
q. Canonical value for the core radius, rc = 5 kpc is used in this 
plot. 

events toward Magellanic Clouds; another, entirely theoret- 
ical, argument is that the cooling times for protogalactic 
halo gas in this region {R < 50 kpc) are an order of mag- 
nitude shorter than the dynamical time (Rees & Ostriker 
1977; White & Rees 1978), thus making a collapse of bary- 
onic structures a likely outcome. We shall extend somewhat 
the discussion of the relevant physics in the Sec. 7. As we 
shall see from the plots in this and subsequent figures, fur- 
ther increase in Rt leads to huge masses of MACHO haloes, 
which are unacceptable from the point of view of BBNS, un- 
less the fraction of galaxies containing MACHO haloes simi- 
lar to the one of the Milky Way is, for some quite mysterious 
reason, very small. 

We also note that MACHOs are incapable of explaining 
galactic dynamics on large scale, and in order to explain dy- 
namical estimates of the Milky Way mass based on satellite 
systems (Kulessa & Lynden-Bell 1992; Zaritsky et al. 1997) 
dark matter in form of either invisible gas or non-baryons 
has to be invoked. As we shall see below, if one chooses to 
accept the BBNS constraints for the baryonic cosmological 
density fls, non-baryonic dark matter in the Milky Way 
halo (and haloes of all normal galaxies) has to be invoked. 
On the other hand, the presence of non-baryonic dark matter 
within MACHO halo itself (i.e. at galactocentric distances 
up to ~ 50 kpc) has the effect of relaxing bounds on the 
MACHO halo, and thus is, in principle, subject to empirical 
verification through detailed comparison of observed opti- 
cal depths with those predicted by the density profile in the 
Eq. (7). 

As far as other components of the total cosmological 
density in baryons Qb are concerned, we adopt the Per- 
sic & Salucci (1992) estimate for the mass of visible baryons 
(i.e. stars, ISM and gas in rich clusters) given by the Eq. (4). 
Further contribution is expected to come from the inter- 
galactic and/or "invisible" galactic gas. At later epochs, 
this is what FHP call "warm gas around galaxies and small 
groups". Since these are presumably the same objects as 
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Figure 11. Masses of MACHO haloes (in units of 10^^ Mq) for 
three different choices of the flattening parameter q, as a function 
of the core radius rc for a fixed truncation radius of Rx = 50 
kpc. Square points represent <? = 1 case, circles q = 0.2 extremely 
flattened case, and crosses the intermediate one {q = 0.6). 



those recently detected as the dominant fraction of the low- 
redshift Lya forest, we shall denote this contribution as flhya 
in further discussion. 

There are basically two ways in which one can discuss 
relationship between Lya clouds and baryonic dark matter: 
(i) direct comparison of their cosmological density in various 
epochs (as in EEC), and (ii) we can consider transformation 
of high-2: Lya clouds into present day MACHOs. We proceed 
with (i), and shall return to the topic (ii) in the next section. 

In order to translate individual galactic masses, as in the 
Eq. (9) into global cosmological density parameter fiiviACHO , 
it is necessary to perform integration over the luminosity 
function (LF), with some assumptions. Beside universality 
of the luminosity function, we have to assume that there 
is no diffuse, intergalactic population of MACHO-like ob- 
jects. This is not just a formal statement - it puts obvi- 
ous constraints on epoch of formation of such objects and 
their degree of clustering. It is natural to speculate that, due 
to dynamical effects, some MACHOs will be ejected from 
the halo during galactic history, thus creating such an in- 
tergalactic population, with its own particular contribution 
to the value of Qb- In this sense, our present picture is not 
completely self-consistent, since it neglects this intergalac- 
tic population of collapsed objects (expression "MACHO" 
is, obviously, inadequate here). In the course of future work, 
we hope to quantify this assumption in detail and, especially, 
demonstrate implications for high-density regions (e.g. rich 
clusters), where "sharing" of the BDM among galaxies may 
have crucial infiuence upon its evolutionary history, and re- 
sult in observable peculiarities (e.g. White & Fabian 1995). 

Cosmological density parameter in such MACHOs re- 
siding in haloes of typical luminous galaxies can be written 
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as 

nMACHO = — / M{L)ip{L) dL, (10) 

Pcrit J 

i'min 

where pcrit is the critical density of Priedmann-Robertson- 

Walker universes, Lmin and Lmax are the minimal and max- 
imal luminosity of galaxies possessing such MACHO haloes, 
respectively, and ip{L) is the universal galaxy LF (Schechter 
1976; Binggeli, Sandage & Tamman 1988; Willmer 1997). 
We use the LF in Schecther's form 

^(L) = ^, (i-) 'cxp(-i-), (11) 

and alternative mass form can be found in Nulsen & Fabian 
(1997). LF parameters are chosen to be (Willmer 1997) 

93* = 2.5 X 10"^ Mpc"^, (12) 

and 

7 = 1.27. (13) 

The influence of their variation on our results is discussed in 
detail below. The standard Schechter luminosity is chosen 
to be L, = 1.0 X 10^° e^°-^^ Lq, i.e. corresponding to 
the absolute B-band magnitude of M, = —19.2 (Willmer 
1997). The Hubble parameter h is chosen to be /i = 0.5 in 
all calculations, except where otherwise mentioned. 

We note from the Eqs. (9) and (11) that a unique value 
of fiMACHO corresponds to each pair of values (g, 7?t). The 
distribution of possible values of this cosmological density 
parameter is shown by the 3-D plot in Fig. 12 for rc = 5 kpc 
and h = 0.5. 

Wo have used fiducial value for the core radius Vc — 5 
kpc. Variation of this quantity in the usual range 5—8 kpc 
(Binney & Tremaine 1987) causes changes in our results of 
6M/M = 6Q/ii < 12 per cent (sec Fig. 11). In addition, we 
have also investigated somewhat unorthodox value rc = 20 
kpc (but, see Gerhard 1999). This is motivated by some re- 
cent indications that the Milky Way rotation curve may be 
satisfactorily explained by nearly homogeneous dark mat- 
ter distribution within a few solar circles (Ninkovic, private 
communication); sec also Fricman & Scoccimarro (1994). 
Also, such a large core-radius would be completely in ac- 
cord with suggestion, originating with N-body simulations 
(Colo & Lacoy 1996), that the density profile becomes signif- 
icantly fiatter then the isothermal one in inner halo regions 
(becoming simultaneously steeper than r^'^ in outermost re- 
gions). The influence of varying core radius on the mass of 
a fiducial halo with our model profile is shown in Fig. 11. 

Variation of other parameters also does not remedy the 
high value of the total mass in MACHOs. For example, Voo 
is only bound from below, by the lAU value of Galactic rota- 
tion at the solar circle of 220 km which was used in these 
calculations. If, as indicated. Milky Way rotation curve rises 
all the way to ~ 3_Ro (Ninkovic, private communication), 
Voo can be as high as 280 km s^^ (Fricman & Scoccimarro 
1994), and the mass M{q, Rt) would be increased for a fac- 
tor 1.62 with corresponding increase in JImacho, which, 
taking into account the bounds in the Eq. (1), is not insignif- 
icant. 

In Figures 12 and 13, we have shown the total MACHO 
+ visible cosmological density vs. the constraints emerging 




Figure 12. The cosmological density parameter Q in MACHOs 
plus visible baryons as a function of the truncation radius, Rt and 
the flattening parameter, q. Canonical value for the core radius, 
rc = 5 kpc is used in this plot, and h = 0.5. Lower bound from 
the primordial nucleosynthesis is also shown, for comparison. 




Figure 13. The same as in Fig. 12, except that the upper nucle- 
osynthctic bound is shown (for the same choice of parameters).. 
It seems clear that we need a significant gaseous baryonic compo- 
nent in the present day universe, quite in accord with the low-z 
Lya forest observations. 

from the BBNS for h = 0.5. We notice that MACHOs within 

50 kpc arc certainly capable of solving the problem of miss- 
ing baryons resulting from comparison of the Eqs. (2) and 
(4). Again, exceptions are very flattened haloes with q ~ 0.2, 
which make them still less appealing possibilities. On the 
other hand, if higher nucleosynthesis estimates of for ex- 
ample Qb ~ 0.077 (Buries & Tytler 1998), are reconfirmed 
by impending observations, there seems to be no alternatives 
to discarding little flattened values q > 0.8 either. 

One should always keep in mind that there is no phys- 
ical reason for assumption that Rt is close to the canonical 
value of 50 kpc; rather, it is just an empirical convenience, 
at least for the time being. Caution suggests to take these 
values (i.e. Rt and the corresponding masses) as lower lim- 
its only, with consequences that fiattening looks even more 
appealing as a way to reduce SImacho- If MACHO haloes 
extend to anything similar to the extent of dynamical haloes 
inferred, for example, by Zaritsky et al. (1997), then reject- 
ing of anything with q > 0.6 is unavoidable. On the other 
hand, it is possible that our reliance on the Occam's razor is 
misleading, and only some fraction / of the mass distribu- 
tion creating potential responsible for the rotational curve is 
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in form of M ACHOs. The rest 1 — / must then be in the form 
of non-baryonic dark matter, if we wish to remedy the high 
fis problem, which leads to a degeneracy, where flattened 
fuU-MACHO halo may contain the same amount of mass 
as non- flattened realistic halo with / < 1. On the other 
hand, optical depth estimates and ratios discussed in the 
Sec. 4 would still be dependent only on the MACHO frac- 
tion, and therefore optical depths are expected to be reduced 
by the factor / and their ratios to be unaffected. This offers 
a further opportunity for improved microlensing statistics, 
which should be able to easily discriminate between values 
of / close to unity and any other significantly smaller value 
(at present, as visible from the Table 1 and comparison with 
Figs. 4 through 6, we can only claim / > 0.1 with reasonable 
certainty) . 



6 GASEOUS CONTENT OF GALAXIES AND 
FLATTENING 

The discovery that a large fraction of low-rcdshift Lya forest 
is associated with normal luminous galEixies (Spinrad et al. 
1993; Lamzetta et al. 1995; Chen et al. 1998) presents a fur- 
ther difficulty for the total baxyonic census, as recognized by 
FHP. This means that at least some debris from the galaxy 
formation epoch remained in the gaseous state till relatively 
late epochs; the question whether this gas (discovered up to 
huge galactocentric distances, with maximal absorption ra- 
dius for L, galaxies being ~ 300 kpc) was partially recycled 
through some galactic stellar population is unimportant in 
this respect. Possible contribution to the baryonic budget is 
enormous; it is the largest (albeit the most uncertain) en- 
try in the list of FHP. In fact, its magnitude is such that 
the BBNS constraints arc seriously jeopardized by a direct 
extension of column-density statistics to the total mass con- 
tained along all linos of sight, an insight which prompted 
RH95 to suggest a significant flattening of these gaseous 
structures. 

In that work, it is shown that constraints following from 
the general formula for the cosmological density of Lya sys- 
tems 

JVmax 

n^^^ = lf!BEl f x-'Nf{N)dN, (15) 

CPcrit J 

(where is the neutral hydrogen column density spanning 
the interval between A^min and Nmax, x neutral gas fraction, 
and f{N) the neutral hydrogen column density distribution) 
coupled with the BBNS bounds leads to inevitable conclu- 
sion of global flattening, if sizes (or coherence) lengths ob- 
tained from double line-of-sight analyses are taken seriously. 
Typical values of Qhya — 0.04 are obtained for typical sizes 
of 100 kpc, in spherical case, from several of their simple 
models, which is, obviously, quite high. If large coherence 
sizes inferred from double linos of sight (e.g. Dinshaw et 
al. 1995) arc characteristic (and they arc in general accord 
with the huge sizes of galactic gaseous haloes obtained by 
Chen et al. [1998]) the cosmological density is even higher. 
The way out is to assume that axial ratio of these struc- 
tures (without entering the question of their physical ori- 
gin and location) is small, and for the most conservative of 
their models, RH95 obtain q^ya < 0.1. It is interesting that 



they suggest clumping of the neutral content as an alterna- 
tive way of decreasing the total mass, a frequent suggestion 
which has not been fully investigated to date (e.g. Mo & 
Miralda-Escudc 1996). 

Conclusions of RH95 are, it should be reemphasized, 
essentially independent of the true nature and location of 
the Lya forest clouds. They are valid for both inter- and 
intragalactic types of absorbers. But the discovery of large 
population of halo absorbers at z < 1 prompts us to ask 
whether aspect ratio of absorbers in RH95 can, in fact, be 
interpreted as the flattening parameter of gaseous haloes. In 
addition, not only the fact that low- and intermediate- 2; ab- 
sorbing clouds preferentially lie in galactic haloes, without 
noticeable morphological segregation (Yahata et al. 1998), 
but also the fact that the covering factor of such haloes was 
found to be close to unity everywhere within the absorb- 
ing radius (Chen et al. 1998), suggests that gaseous haloes 
should be flattened at late epochs. It should be noted that 
flattening of Lya absorption systems was much earlier pro- 
posed, for different reasons, by Barcons & Fabian (1987) and 
Milgrom (1988). 

Unfortunately, exact knowledge of the baryonic content 
of the Lya absorbing clouds requires certain knowledge on 
their ionization structure, which is still very elusive. The 
value of metagalactic ionizing background, which is the only 
always operating ionizing source, is still painfully uncertain 
even in the local universe and at low redshift, and the more 
so at high-z (Bajtlik, Duncan & Ostriker 1988; Kulkarni & 
FaU 1993; Vogel et al. 1995; Donahue, Aldering & Stocke 
1995). The presence of internal ionizing sources, inferred in 
some local intergalactic clouds (Donahue et al. 1995; Bland- 
Hawthorn et al. 1995) is also quite uncertain. Finally, geo- 
metric properties of ensembles of clouds, i.e. their dumpiness 
and global flattening, are still only speculative. 

The comparison of FFG with the mass estimate of the 
Lya forest of Weinberg et al. (1997) is interesting, espe- 
cially in the view of possible absence of mixing between the 
two types of unseen baryonic matter. Weinberg et al. (1997) 
value, quoted by FFG, 

nLyc = 0.02ft"^ (16) 

is not at all especially high when considered within a "fam- 
ily" of closure fractions obtained for fiLya 

Higher values are required for high-« intergalactic Lya 
forest by many models, e.g. Bi & Davidsen (1997) suggest 
f^Lya = 0.025 hr'^ , adopting the values of metagalactic ion- 
izing flux from Haardt & Madau (1996). Alternatively, one 
may wish to reduce the value of fiLya , but at a price of hav- 
ing significantly different ionizing background (since models 
are able only of constraining choices for O'^/Juv)- Further- 
more, any increase in the cosmological bias would tip the 
scales toward larger contribution of gas in comparison to 
the visible matter assembled in stars and "normal" lumi- 
nous galaxies. 

In our opinion, the association of significant fraction (if 
not all) low-z Lya forest with galactic haloes does give ad- 
ditional credence to FFG conjecture about MACHOs being 
distinct baryon reservoir from Lya forest. At least this is so 
after some particular epoch, which wc shall denote by Zhd, 
when bulk of today's MACHOs was formed out of gas-rich 
protogalactic fragments, and which may be called the epoch 
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of baryonic decoupling. 

Gaseous density profile and global shape should, in prin- 
ciple, follow profile and shape of the underlying dark matter 
distribution, which is supposed to be dynamically dominant 
(Nulsen 1986). This is the necessary link between arguments 
concerning flattening of gaseous structures around galaxies 
and flattening of MACHO halo. Of course, exact amount of 
flattening (as well as the exact density profile) would be dif- 
ferent because dark matter is, in contradistinction to gaseous 
matter, assumed to be dissipationless (no matter whether it 
is mainly MACHO or non-baryonic elementary particle), but 
the difference is not expected to be very large on the scales 
of ~ 50 kpc. If there is a significant infall from the halo to 
the disc, the constraints could be much tighter. In this re- 
spect, it is important to mention that an argument of this 
type was advanced by Sancisi & van Albada (1987), namely, 
that recent accretion of gas found in the optical plane at 
outer fringes of spiral galaxies already implies flattening of 
the dynamical halo in these regions. 

Conclusions of FFG actually receive multifold support 
from considerations of low-redshift Lya forest: 

1. Probable overestimate of the f^Lya in N-body simula- 
tions and other models of unclustered population of Lya 
clouds makes it easier to accommodate large Omacho ■ 

2. The very fact that we perceived much stronger Lya ab- 
sorption at recent epochs than expected on the basis of 
naive extrapolation supports the conclusion that Lya 
clouds arc rather decoupled from star-formation histo- 
ries (at least after some fiducial epoch za), so the dis- 
cussion in the Sec. 3 of FFG is justified. 

In addition, recent observational indications that Milky 
Way still possesses extended gaseous halo with densities 
■~ 10^'' cm^"^ at galactocentric distances ~ 50 kpc (i.e. simi- 
lar to those discussed with respect to the MACHO halo) un- 
derlie the necessity of having large fraction of dark baryons 
in gaseous form at present day (Weiner & Williams 1996). 
This further narrows available range for Omacho, and sug- 
gests that all possibilities to reduce it should be explored, 
flattening being the simplest one of them and most in line 
with the principle of economy of hypotheses. 

Various baryonic components are represented in the 
Q, — h diagrams in Figs. 14 and 15. If the Lya mass es- 
timated by Weinberg et al. (1997) is correct, and the as- 
sumption of FFG about essential decoupling of the baryonic 
contents of Lya forest and MACHOs after some initial pe- 
riod, high values of h > 0.8 seem to be highly implausible for 
both flattened and unflattened MACHO haloes. For spher- 
ical haloes, we ran into troubles for almost all values of h. 
It is marginally acceptable for h = 0.5, but it is inconsistent 
with any higher values (again, wo should keep in mind that 
there is no physical reason for the truncation of MACHO 
halo at the LMC distance, and many arguments that dy- 
namical haloes extend much further). Flattened haloes, on 
the other hand, are quite securely within the BBNS margin 
for h ~ 0.5. 

7 DISCUSSION 

Mass considerations could, in principle, lead us to an impor- 
tant clue in solving the puzzle of the fate of the halo gas. 
All scenarios of hierarchical structure formation (e.g. White 
& Rees 1978; Navarro & White 1994) have massive dark 
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Figure 14. Various components of the baryonic cosmological 
density (visible, Lya absorbing gas and MACHOs) compared with 
the BBNS bounds for varying Hubble parameter h. In this Figure, 
the minimal mass of the MACHO halo within 50 kpc in our dis- 
cussion, corresponding to q = 0.2, is talcen. The sum of all three 
main baryonic components is shown as the dash-double-dotted 
line labelled with S^^™^"*. 

haloes in place by z ~ 2 — 3, which form gravitational poten- 
tial wells accreting diffuse baryonic matter (Mo & Miralda- 
Escude 1996). It is almost certain that at high and inter- 
mediate redshift exactly this gas dominates cosmic baxyonic 
budget (RH95; FHP). Its subsequent history is still shrouded 
in mystery, but the fact that estimates of the cosmological 
density fraction in MACHOs in both FFG and present work 
point out that at some point in the galactic history the tran- 
sition between the halo gas (or its significant fraction) and 
collapsed objects, like the present day MACHOs occurred. 

As the phase transition of baryonic dark matter we 
denote the process of transformation of pregalactic/proto- 
galactic gas into MACHOs. It seems obvious that transition 
from the diffuse (gas) to the collapsed (MACHO) phase of 
the EDM must have occurred at some point in the history of 
the universe. This process, whatever form and in whichever 
epoch it took place, is of the crucial importance for un- 
derstanding the evolution of baryonic content of the uni- 
verse. Parenthetically, this is return to the authentic physical 
meaning of the concept of phase of matter, since Lya clouds 
and MACHOs do possess different symmetry properties. 

This transition from gaseous to collapsed phase can pro- 
ceed in several ways. While exact behavior is still too diffi- 
cult to investigate in detail, since it depends on the initial 
dynamical and chemical conditions, we can sketch several 
possibilities. 
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Figure 15. The same as in Fig. 14, for the maximized MACHO 
cosmological density (non-flattened g = 1 case). 

An early transition probably occurs through Popula- 
tion III stars, and in this scenario, MACHOs are stellar 
remnants (Carr 1994; FFG). Only detailed chemical mod- 
elling cam show whether this is viable. Apart from possible 
overproduction of metals, a further problem with this pic- 
ture is the necessity for the Pop III initial mass function 
(IMF) to be significantly different from that observed today 
in the ISM. This is not only to avoid conflicts with results of 
deep searches (Richstone et al. 1992; Flynn, Gould & Bah- 
call 1996), but also to avoid problems with massive stellar 
black holes (disruption of stellar discs, excessive X-ray emis- 
sion due to accretion, etc.). 

Otherwise, it can proceed via cold gas, passing through 
a phase similar to that depicted in models of Pfenniger, 
Combes & Martinet (1994), Pfenniger & Combes (1994), 
Gerhard & Silk (1996) or Walker (1998) and Walker & War- 
die (1998). It is not clear how the IMF in such cases can be 
restricted to low-luminosity objects. On the other, it may 
be significant that Pfenninger et al. (1994) models do pre- 
dict extreme flattcrnng of baryonic haloes (or thick discs), 
composed of fractal distribution of small, sub-Jeans mass 
cloudlets (Pfenniger & Combes 1994). 

Finally, gas may form MACHOs in an early cooling flow, 
similar to those we perceive today in rich clusters and around 
isolated giant elliptical galaxies (Fabian, Nulsen & Canizares 
1982; Sarazin 1988; Fabian et al. 1986; Nulsen, Johnstone & 
Fabian 1987). This is scenario of Nulsen & Fabian (1995; 
1997). It has the advantage that it creates conditions favor- 
able to low-luminosity objects in a natural way (Sarazin & 
O'Connell 1983; Sarazin 1986; Ferland, Fabian & Johnstone 
1994; Fabian & Nulsen 1994). Further virtue of this picture 



is that it presents continuity with cooling of hot, virial haloes 
recently invoked in order to explain metal and Lyman-limit 
absorption systems (Mo 1994; Mo & Miralda-Escude 1996). 

As Nulsen & Fabian (1997) point out, cooling flow sce- 
nario seems to be capable of solving the fundamental puzzle 
of origin of morphological difference between galaxies. This 
is even more important in view of very recent QSO absorp- 
tion results indicating that there is no signiflcant difference 
between absorption properties of low- and intormediate- 
redshift early- and late- type galaxies (Yahata et al. 1998). 
Since Lya absorption naturally samples larger galactocen- 
tric distances than optical surveys which produced all ex- 
isting morphological classiflcations, it follows that outer re- 
gions, where the dynamical timescalc is longer, are similar 
for all haloes of comparable total (i.e. baryonic -|- nonbary- 
onic) mass. Since thermal instabilities, cooling and infall are 
necessary consequences of an early virialization (e.g. White 
& Rees 1978), one does flnd present-day cooling flows with 
their large mass deposition, a viable model for the evolu- 
tion of gaseous haloes. A parcel of infalling cooled gas, as 
pointed out by Binney (1995) and Mo & Miralda-Escude 

(1996) , can have only two fates: it may fall into the disc and 
be supported by angular momentum, effectively becoming 
part of the galactic ISM, or it can experience sufficient cool- 
ing before reaching the disc, or the centre of the halo, and - 
possibly passing through transient molecular phase - create 
stars and/or brown dwarves. Those two paths are likely to 
produce entirely different global EDM structures; the former 
would create (and help to sustain) spiral galaxies, and the 
latter would apply to triaxial systems. The detailed models 
along these lines, extending the view of Nulsen & Fabian 

(1997) , may present a significant stop forward in explaining 
the morphological differences between galaxies. 

In any case, it is our hope that future detailed mod- 
elling, coupled with observational breakthroughs in both 
microlensing and halo absorption systems will answer the 
question on the empirical value of the baryonic fraction of 
gaseous haloes, which can be written as 

„ _ MmACHO + Mgas + Myis 

- m;;, • ^^^> 

This parameter is crucial for the theories describing post- 
virialization cooling and early infall of gas. Mo & Miralda- 
Escude (1996) choose fg = 0.05, value that is in general 
agreement with orthodox assumptions of the BBNS, but oth- 
erwise remains a free parameter of the model. Large cosmo- 
logical constant, for example, like recently popular models 
with Qa ~ 0.6 (e.g. Schmidt et al. 1998), although obviously 
not influencing the cooling of protogalactic gas, will manifest 
itself indirectly through a significant increase in fg. 

Pietz et al. (1998) show that X-ray emission from the 
inner Milky Way gaseous halo is best reproduced in a flat- 
tened model with q ~ 0.3, although other flattened models 
give better agreement with the data than the spherical mod- 
els. 

The flattened shape of the dark halo is not orfly im- 
portant in connection with the EDM, but also with non- 
baryonic dark matter such as massive neutrinos. Decaying 
dark matter theory (Sciama 1993), which is based on the 
assumption that the major constituent of the dark mass in 
haloes of spiral galaxies is a massive decaying neutrino, with 
the mass of ~ 30 eV and lifetime ~ 2 x 10^^ s, requires that 
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the haloes are extremely flattened - the halo of the Galaxy 
is flattened so the axial ratio is g = 0.2 (Sciama 1990, 1997). 
This is a consequence of the eSort to reconcile the observed 
electron densities obtained from the pulsar dispersion mea- 
sure data (Nordgren, Cordes & Terzian 1992) according to 
which Tie ~ 0.03 cm~^. DDM theory would give the value of 
He ~ 0.017 under the assumption of spherical halo. Agree- 
ment is thus obtained with significant flattening of the halo 
that reduces the scale height by the factor ~ 4 to ~ 2 kpc. 
One can notice that this flattening could be achieved by 
extending the mass models of the Galaxy by Dehnen & Bin- 
ney (1998). However, recently, an attempt has been made 
to show that in the case of the Galaxy the axis ratio q is 
q — 0.75 ± 0.25 thus ruling out cold molecular gas and de- 
caying massive neutrino as viable dark matter candidates 
(Oiling & Merrifield 1998). It is somewhat beyond the scope 
of the present paper, and we only notice that the aforemen- 
tioned result with this, rather high, value of the parameter 
q is attained if the galactocontric distance is = 7± 1 kpc, 
which is a rather unorthodox value. On the other hand, mi- 
crolensing optical depths indicate that "standard" versions 
of these theories, requiring q ~ 0.2 are no longer viable. This 
conclusion can be avoided if MACHOs are dynamically in- 
significant within 50 kpc. However, it can be shown that 
one can easily accommodate a much larger value of q, i.e. 
q ~ 0.6, into the DDM theory without significant changes of 
the theory's fundamental parameters such as mass and life- 
time of the decaying neutrino (Samurovic & Cirkovic 1999). 



8 CONCLUSIONS 

On the basis of still scarce empirical data, and still undevel- 
oped and unsophisticated theoretical models, one can, how- 
ever, draw some important inferences regarding the global 
shape of haloes of spiral galax:ies, taking the Milky Way as 
a prototype L ~ L, galaxy at zero-redshift. It seems that a 
whole array of different arguments point to an oblate grav- 
itating dark halo, which, in turn, causes flattening of other 
types of haloes, like stellar and gaseous (Lya-absorbing) 
haloes. Our conclusions are thus summarized as follows: 

1. The set of measured optical depths for microlensing, 
although still statistically incomplete, strongly indicates 
moderately flattened haloes. 

2. We reconfirm the conclusion of FFG that Omacho is 
very high in all plausible cases, specifically, 

f^MACHo/f^s =0.1 — 1. 

For spherical or little-fiattened haloes, it is, in fact, un- 
pleasantly close to Jls, as obtained from the primordial 
nucleosynthesis and flattening seems to bo the simplest 
remedy for this situation. 

3. Review of relevant literature reveals a multitude of ar- 
guments for flattened haloes of the Milky Way and other 
spiral galaxies. We find the arguments based on the to- 
tal mass of Lya forest clouds especially convincing, in 
conjunction with strong arguments for association of 
low and intermediate-redshift Lya forest with normal 
galaxies. 

4. We find the FFG conclusion that baryons in MACHOs 
and Lya forest are essentially decoupled is strengthened 
on the basis of low- 2; absorption studies and indications 



of extended gas around present-day galaxies. 
5. Baryonic census clearly favors low values of the Hubble 

constant, essentially irrespectively of flattening. 
Further investigations, especially of the epoch of MACHO 
formation, probed by early damped Lya and similar gas- 
rich systems will be necessary to completely clear the pic- 
ture of ramification of the baryonic matter in the universe 
into diffuse and collapsed components. This, coupled with 
microlensing advances, should bo able to finally solve the 
problem of baryonic component of dark matter and unify 
several branches of astrophysical research into a coherent 
picture of the evolution of matter, as we know it, in the 
universe. 



ACKNOWLEDGEMENTS 

The authors are happy to express their gratitude to Dr. Geza 

Gyuk for stimulating comments and encouragement. Useful 
discussions with Dr. Slobodan Niukovic and helpful remarks 
by Dr. Giuliano Giuricin are also acknowledged. SS acknowl- 
edges the financial support of the University of Trieste where 
part of this work was carried out. 



REFERENCES 

Alcock, C. et al., 1996, ApJ, 461, 84 
Alcock, C. et al., 1997a, ApJ, 479, 119 
Alcock, C. et al., 1997b, ApJ, 486, 697 
Alcock, C. et al., 1997c, ApJ, 491, Lll 
Ansari, R. et al., 1996, A & A, 314, 94 
Aubourg, E. et al., 1993, Nat, 365, 623 
Bahcall, J.N., 1984, ApJ, 287, 926 
Bahcall, J. N. 1986, ARA&A, 24, 577 
Bahcall, J.N., Spitzcr, L., 1969, ApJ, 156, L63 
Bajtlik, S. Duncan, R. C, Ostrikcr, J. P., 1988, ApJ, 327, 570 
Barcons, X., Fabian, A. C, 1987, MNRAS, 224, 675 
Bergeron, J., Boisse, P. 1991, A& A, 243, 344 
Bi, H., Davidsen, A. F., 1997, ApJ, 479, 523 
Binggeli, B., Sandage, A., Tammann, G. A., 1988, ARA&; A, 26, 
509 

Blnney, J., 1995, Oxford Univ. preprint OUTP/95/09A 
Binncy, J., Gcrliard, O.E., Spcrgcl, D., 1997, MNRAS, 288, 365 
Binney, J., Gerhard, O.E., Stark, A. A., Bally, J., Uchida, K.I., 

1991, MNRAS, 252, 210 
Binney, J., May, A., Ostriker, J.P., 1987, MNRAS, 226, 149 
Binney, J., Merrifield, M., 1998, Galactic Astronomy, Princeton 

Univ. Press, Princeton, NJ 
Binney, J., Tremaine, S., 1987, Galactic Dynamics, Princeton 

Univ. Press, Princeton, NJ 
Bland-Hawthorn, J., Ekers, R. D., van Bruegel, W., Koekemoer, 

A., Taylor, K., 1995, ApJ, 442, L77 
Bristow, P. D., Phillipps, S., 1994, MNRAS, 267, 13 
Buries, S., Tytler, D., 1998, ApJ, 499, 699 
Canuto, v., 1978, MNRAS, 184, 721 
Carr, B. J., 1994, ARA& A, 32, 531 

Chen, H-W., Lanzetta, K.M., Webb, J.K., Barcons, X., 1998, ApJ, 
498, 77 

Cole, S.. Lacey, C., 1996, MNRAS, 281, 716 

Combes, F., Boisse, P., Mazure, A., Blanchard, A., 1995, Galaxies 

and Cosmology, Springer- Verlag, Berlin 
Crotts, A. P. S., 1992, ApJ, 399, L43 
Crotts, A. P. S.. Tomaney, A. B., 1996, ApJ, 473, L87. 
De Riijula, A., Jetzer, Ph., Masso, E., 1992, A& A, 254, 99 
Dehnen, W., Binney, J., 1998, MNRAS, 294, 429 



Flattened Galactic Haloes and Baryonic Dark Matter 13 



Derue. F. ct a l. (EROS collaboration), 1999, preprint astro- 



ph/9903209 



Dinshaw, N., Impey, C. D., Foltz, C. B., Weymann, R. J., Morris, 

S. L., 1995, Nat, 373, 223 
Eichler, D., 1996, ApJ, 468, 75 

Fabian, A. C, Arnaud, K. A., Nulsen, P. E. J., Mushotzky, R. F., 

1986, ApJ, 305, 9 
Fabian, A. C, Nulsen, P. E. J., 1994, MNRAS, 269, L33 
Feast, M., Whitelock, P., 1997, MNRAS, 291, 683 
Ferland, G. J., Fabian, A. C, Johnstone, R. M., 1994, MNRAS, 



m, m 

•Fiylds, D. D., FmusH, K., GialT, D. 3., 1998, NewA, 3, 347 (FFG) 
Flynn, C, Gould, A., Bahcall, J. N., 1996, ApJ, 466, L55 
Freudenreich, H.T., 1998, ApJ, 492, 495 
Frieman, J., Scoccimarro, R., 1994, ApJ, 431, L23 
Fukugita, M., Hogan, C.J., Peebles, P.J.E., 1998, ApJ, 503, 518 
(FHP) 

Gates, E., Gyuk, G., Turner M. S.. 1995. ApJ . 449, 123 
Gerhard, O., 1999, preprint ptro-ph/9902247 
Gerhard, O., Silk, J., 1996, ApJ, 472, 34 

Gilmore, G., Wyse, R. F. G., Kuijken, K., 1989, ARA& A, 27, 
555 

Gnedin, N. Y., Ostriker, J. P., 1992, 400, 1 
Gould, A., 1994, ApJ, 435, 573 
Gould, A., 1996, PASP, 108, 465 

Gould, A., Bahcall, J.N., Flynn, C., 1997, ApJ, 482, 913 

Graff, D. S., Freese, K., 1996, ApJ, 467, L65 

Griest, K., et al., 1991, ApJ, 372, L79 

Gyuk, G., Holder, G.P., 1998, MNRAS, 297, L44 

Hata, N., Scherrer, R. J., Steigman, G., Thomas, D., Walker, T. 

P., 1996, ApJ, 458, 637 
Hegyi, D. J., Olive, K. A., 1986, ApJ, 303, 561 
Henriksen, R. N., Widrow, L. M., 1995, ApJ, 441, 70 
Hu, E. M., Huang, J.-S., Gilmore, G., Cowie, L. L., 1994, Nat, 

371, 493 

Jezter, P., Striissle, M., Wandeler, U., 1998, A & A, 336, 411 
Kiraga, M., Paczyhski, B., 1994, ApJ, 430, LlOl 
Kulessa, A., Lynden-Bell, D., 1992, MNRAS, 255, 105 
Kolb, E.W., Tkachev, 1995, Phys. Rev. D, 50, 769 
Koopmans, L. V. E., de Bruyn, A. G., Jackson, N., 1998, MNRAS, 
295, 534 

Kovalevsky, J., 1998, ARA& A, 36, 99 
Kuijken, K.. Gilmore, G., 1991, ApJ, 367, L9 
Lanzetta, K. M., Bowen, D., Tytler, D., Webb, J. K., 1995, ApJ, 
442, 538 

Lee, M.G., Freedman, W., Mateo, M., Thompson, I., Roth, N., 

Ruiz, M.-T., 1993, AJ, 106, 1420 
Mc Uier. Y.. Bernarde au. F., Van Waerbeke, L., 1998, preprint 



istro-ph/ 9802005 



Merrificld, M.M., 1992, AJ, 103, 1552 
Milgrom, M. 1988, A & A, 202, L9 
Mo, H. J., Miralda-Escude, J., 1996, ApJ, 469, 589 
Navarro, J. F., White, S. D. M., 1994, MNRAS, 267, 401 
Nakamura, T., Kan-ya, Y., Nishi, R., 1996, ApJ, 473, L99 
Ninkovic, S., 1985, Ap& SS, 110, 379 
Ninkovic, S., 1991, Bull. Obs. Astron. Belgrade, 143, 49 
Ninkovic, S., Petrovskaya, I. V., 1992, AZh, 69, 926 
Nordgren, T. E., Cordes, J. M., Terzian, Y., 1992, AJ, 104, 1465 
J., 1986, MNRAS, 221, 377 

J., Johnstone, R. M., Fabian, A. C., 1987, PASA, 



Nulsen, P. E. 
Nulsen, P. E 

7, 132 
Nulsen, P. E. 
Nulsen, P. E. 



J., Fabian, A. C., 1995, MNRAS, 277, 561 
J., Fabian, A. C., 1997, MNRAS, 291, 425 
Oiling. R.P., 1995, PhD thesis, Columbia Univ. 
Oiling, R.P., Merrifield, M.R., 1998, in Zaritsky D., ed.. Galactic 
Halos: A UC Santa Cruz Workshop, ASP Conf. Ser. Vol. 136, 
Astron. Soc. Pac, San Francisco, p. 219 
Oiling, R.P., Merrifield, M.R., 1998, MNRAS 297, 943 



Oort, J. H., Plant, L., 1975, A & A, 41, 71 
Paczyhski B., 1986, ApJ, 304, 1 

Paczyhski B. et al. (OGLE Collaboration), 1994, ApJ, 435, Lll 
Paczyhski B., 1994, Acta Astron. 44, 235 

Palanque-Delabrouille, N. et al. (EROS Collaboration), 1998, A 
& A, 332, 1 

Peebles, P. J.E., 1993, Principles of Physical Cosmology, Princeton 

Univ. Press, Princeton, NJ 
Persic, M., Salucci, P., 1992, MNRAS, 258, 14p 
Persic. M.. Salu cci. P., 1998, MNRAS, in press (preprint a,stro- 



ph/9806215) 



Pfenniger, D., Combes, F., Martinet, L., 1994, A & A, 285, 79 

Pfenniger, D., Combes, F., 1994, A & A, 285, 94 

Pietz, J., Kerp, J., Kalberla, P. M. W., Burton, W. B., Hartmann, 

D., Mebold, U., 1998, A & A, 332, 55 
Pitts, E., Tayler, R. J., 1997, MNRAS, 288, 457 
Ranch, M., Haehnelt, M. G., 1995, MNRAS, 275, L76 (RH95) 
Rees, M. J., Ostriker, J. P., 1977, MNRAS, 179, 541 
Renault, C. et al. (EROS Collaboration), 1997, A & A, 324, 69 
Rhoads, J.E., Malhotra, S., 1998, ApJ, 495, L55 
Richstone, D., Gould, A., Guhathakurta, P., Flynn, C, 1992, 

ApJ, 388, 354 
Roulet, E., MoUerach, S., 1997, Phys. Rep., 279, 68 
Sackett, P. D., Sparke, L. S., 1990, ApJ, 361, 408 
Sackett, P. D., Gould, A., 1993, ApJ, 419, 648 
Sackett, P. D., Rix, H.-W., Jarvis, B. J., Freeman, K. C, 1994, 

ApJ, 436, 629 
Sackett, P. D., 1997, ApJ, 483, 103 

Samurovic, S., Cirkovic, M.M., 1999, Journal of Research in 

Physics, in press 
Sancisi, R., van Albada, T. S., 1987, in Kormendy J., Knapp 

G.R., eds., Proc. lAU Symp. 117, Dark Matter in the Universe, 

Reidel, Dordrecht, p. 67. 
Sarazin, C. L., 1988, X-Ray Emission from Clusters of Galaxies, 

Cambridge Univ. Press, Cambridge 
Sarazin, C, O'Connell, R. W., 1983, ApJ, 268, 552 
Schechter, P., 1976, ApJ, 203, 297 
Schmidt, B.P. et al., 1998, ApJ, 507, 46 
Sciama, D. W., 1990, MNRAS, 244, Ip 

Sciama, D.W., 1993, Modern Cosmology and the Dark Matter 
Problem, Cambridge Univ Press. Cambridg e 



astro-ph/9704081 



Sciama, D.W., 1997, preprint 
Spinrad, H. et al., 1993, AJ, 106, 1 
Stanek, K.Z., 1995, ApJ, 441, L29 

Steigman, G., Tkachev, I., 1 998. preprint astro-ph/980300j 

Turner, M.S., 1996, preprint :istro-ph/9610158 

Walker, M. A., 1998, preprint ptro-ph/9807236 

Walker, M., Wardle, M., 1998, ApJ, 498, L125 

Walker, T. P., Steigman, G., Schramm, D. N., Olive, K. A., Kang, 

H.-S., 1991, ApJ, 376, 51 
Wasserman, I., Salpeter, E. E., 1994, ApJ, 433, 670 
Weinberg, D. H., Miralda-Escude, J., Hernquist, L., Katz, N., 

1997, ApJ, 490, 564 
Weiner, B. J., Williams, T. B., 1996, AJ, 111, 1156 
White, S. D. M., Rees, M. J., 1978, MNRAS, 183, 341 
White, D. A., Fabian, A.. C, 1995, MNRAS, 273, 72 
White, M., Viana, P. T. P., Liddle, A. R., Scott, D., 1996, MN- 
RAS, 283, 107 
Willmer, C. N. A., 1997, AJ, 114, 898 

Wyse, R. F. G., Gilmore, G., 1989, Comments Astrophys. 13, 135 
Yahata, N., Lanzetta, K. M., Web b. J. K.. Barcons. X., 1998, 



Blois conf. proceedings, preprint astro-ph/980929(: 



Yang, J., Turner, M. S., Steigman, G., Schramm, D. N., Olive, K. 

A., 1984, ApJ, 281, 493 
Zaritsky, D., Smith, R., Frenk, C, White, S. D. M., 1997, ApJ, 

478, 39 



14 S. Samurovic, M. M. Cirkovic and V. Milosevic-Zdjelar 

This paper has been produced using the Royal Astronomical 
Society /Blackwell Science T^^X macros. 



